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FOREWORD 


This  report  presents  essentially  the  sane  material  as  the 
Master's  Thesis,  "An  Experimental  Investigation  of  the  Effect  cf 
Electrically  Induced  Controlled  Frequency  Perturbations  on 
Boundary-Layer  Transition,"  by  Jeffrey  Ketcham.  The  thesis  was 
presented  in  partial  fulfillment  of  the  requirements  for  the 
Master  of  Science  in  Mechanical  Engineering  at  The  Ohio  State 
University. 
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Office  -  Durham.  At  Battelle  Institute  specific  thanks  are  given 
to  J.  E.  Gates,  W.  E.  Kortear,  R.  Badertscher,  V.  0.  Hoehne,  and 
D.  Aithen,  for  their  assistance  in  bringing  this  project  to  fruition. 
The  junior  author  also  wishes  to  express  his  thanks  to  Dr.  Henry  R. 
Velkoff  for  his  suggestion  of  topic,  persistent  encouragement  and 
help  in  the  conduct  of  this  research. 
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ABSTRACT 

An  experimental  investigation  was  conducted  to  determine  the 
effects  of  corona  discharge  on  boundary -layer  transition.  A  flat  plate 
was  mounted  in  Battelle  Memorial  Institute's  low-turbulence  wind  tunnel. 
The  location  of  transition  was  determined  using  a  near-surface  impact 
pressure  probe.  Boundary-layer  excitation  was  accomplished  by  suspend¬ 
ing  a  series  of  corona  wires  parallel  to  the  plate  outside  the  boundary 
layer.  It  was  found  that,  under  certain  conditions,  the  transition 
point  could  be  moved  downstream  upon  application  of  a  steady  potential 
to  the  wires.  Superposition  of  a  pulsating  potential  did  not  influence 
the  results  over  a  wide  range  of  pulsation  frequency. 
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NOMENCLATURE 


Symbol 

Units 

A 

Amplitude  of  boundary-layer  oscillation 

ft 

c 

r 

Velocity  of  boundary-layer  oscillation 

ft/sec 

E 

Electric  field  strength 

volts/ft 

£ 

Frequency 

cps 

K 

Ion  mobility 

2 

ft  /volt-sec 

P 

Total  pressure 

lb/ft2 

P' 

Fluctuating  component  of  pressure 

lb/ft2 

P 

Steady  component  of  pressure 

lb/ft2 

R 

Gas  constant 

ft-lb./lb  -°R 
f  m 

Re 

Reynolds  Number 

dimensionless 

t 

Time 

sec 

u 

Total  velocity  in  x  direction 

ft/sec 

u' 

Fluctuating  velocity  in  x  direction 

ft/sec 

U 

Steady  velocity  in  x  direction 

ft /sec 

U 

m 

Maximum  velocity  in  x  direction 

ft/sec 

V 

Total  velocity  in  y  direction 

ft/sec 

v' 

Fluctuating  velocity  in  y  direction 

ft/sec 

V 

Steady  velocity  in  y  direction 

ft/sec 

\ 

Total  ion  velocity 

ft/sec 

V 

a 

Air  velocity  in  corona  wind 

ft/sec 

w 

Total  velocity  in  z  direction 

ft/sec 

w' 

Fluctuating  velocity  in  z  direction 

f t/sec 
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if 
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i 

i 

i 

i 
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NOMENCLATURE  (Continued) 

Symbol 

Units 

N 

Steady  velocity  in  z  direction 

ft/sec 

X 

Distance  along  plate 

ft 

y 

Distance  normal  to  plate 

ft 

z 

Distance  normal  to  x  and  y 

ft 

or 

Reciprocal  of  oscillation  wavelength 

a-1 

»i 

Amplification  factor  for  oscillations 

-1 

sec 

»r 

Circular  frequency  of  oscillations 

rad/scc 

6 

Boundary-layer  thickness 

ft 

* 

6 

Momentum  thickness 

ft 

T| 

Similarity  parameter 

dimensionless 

X 

Oscillation  wavelength 

ft 

Viscosity 

lb  /ft-sec 
m 

V 

Dynamic  viscosity 

2 

ft  /sec 

P 

Density 

lb  /ft3 
m 

♦ 

Stream  function 

ft^/sec 

Subscripts 

L 

Laminar 

T 

Turbulent 

tr 

Transition 
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INTRODUCTION 

Investigations  of  the  transition  of  laminar-to-turbulent 

1* 

flow  date  back  to  the  pipe  and  channel  studies  of  0.  Reynolds  in 
1883.  Such  work  was  prompted  by  the  increased  pressure  gradient 
required  to  maintain  flow  in  a  turbulent  regime.  No  less  significant 
today  are  the  effects  of  boundary-layer  transition  on  the  performance 
of  airfoils,  heat  exchangers,  and  essentially  all  devices  concerned 
with  fluid  motion,  for  associated  with  transition,  are  marked  changes 
in  surface  friction  drag  and  heat-transfer  rate.  While  many  investi¬ 
gations  over  the  years  have  contributed  to  the  characterization  of 
turbulent  flow,  much  is  yet  to  be  learned  regarding  the  detailed  phys¬ 
ical  mechanisms  of  transition  and  possible  means  for  its  control. 

It  is  well  known  that  the  transition  phenomenon  is  sensitive 

to  a  number  of  flow  conditions.  Some  of  these  are  given  below  in 
2 

Table  1  .  Additionally,  it  has  been  experimentally  demonstrated  that 
transition  may  be  forestalled  by  employing  a  flexible  surface  skin,  or 
initiated,  by  a  tripping  device  or  sonic  excitation.  Flexible  skin  for 
delaying  transition  and  reducing  surface  drag  was  first  proposed  by 
Kramer  as  a  result  of  studies  on  the  hydrodynamic  performance  of 

3 

dolphins.  Kramer's  experimental  results  with  flexible  surface  coat¬ 
ings  in  water  showed  that  reduction  of  skin  friction  coefficients  by 
a  factor  of  two  is  possible  under  certain  conditions.  Prandtl  first 
demonstrated  the  effect  of  a  tripping  wire  in  artificially  inducing 


*  Superscripts  denote  items  appearing  in  the  List  of  References. 


2 


transition  on  a  sphere  in  1914.  More  recent  investigations  have  illu- 
strated  the  transitional  effects  of  Jet  engine  noise  ,  loudspeakers  in 
the  flow^,  and  sonic  excitation  at  the  fluid-solid  interface**.  An 
important  consideration  when  predicting  transition  is  that  superposi¬ 
tion  of  the  factors  which  influence  it  is  not  possible,  and  in  each 
flow  case,  the  dominant  effects  must  be  isolated. 

TABLE  1.  FACTORS  INFLUENCING  TRANSITION 


Factor 

Transition  Promoted  By 

Free  stream  turbulence 

Increasing 

Pressure  gradient 

Positive 

Roughness 

Increasing 

Heat  transfer 

Heating 

Surface  curvature 

Increasing  (concave) 

The  above  discussion  points  to  the  fact  that  the  mechanisms 
governing  transition  are  indeed  delicate  and  that  the  external  energy 
involved  to  promote  or  postpone  transition  is  small  compared  to  that 
required  to  maintain  the  basic  flow.  This  observation  in  conjunction 
with  the  theory  of  laminar-boundary-layer  oscillations  forms  the  basis 
for  the  following  thesis.  Boundary-layer  transition  may  possibly  be 
influenced  by  an  electrical  perturbation  of  the  proper  intensity,  fre¬ 
quency,  and  wavelength.  This  thesis  is  developed  in  relation  to  the 
effects  of  corona  discharge  on  the  boundary-layer  transition  of  steady- 
state  incompressible  air  flow  over  a  flat  plate. 
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THEORY  OF  BOUNDARY-LAYER  OSCILLATIONS 


Boundary-layer  stability  theory,  as  ultimately  developed  by 
Tollmein  and  Schlichting,  is  based  on  a  method  of  small  disturbances 
which  is  characterized  by  the  following  premises. 

(1)  Laminar  flow  is  affected  by  small 
disturbances. 

(2)  The  disturbances  are  initially  random  in 
nature  and  are  inherent  in  the  flow,  by 
virtue  of  free  stream  turbulence,  or 
boundary  surface  conditions. 

(3)  If  the  disturbances  decrease  with  time, 
the  flow  will  remain  stable  and  laminar. 

(A)  If  the  disturbances  grow  with  time,  the 

flow  will  become  unstable  with  the  possi¬ 
bility  of  transition  to  turbulence. 

(5)  A  critical  Reynolds  Number  exists  below 
which  no  disturbances  will  be  amplified. 

The  main  lines  of  the  theory  of  stability  are  outlined  in  Schlichting's 
text^  and  discussed  by  Schubauer  and  Skramstad  in  Reference  5.  The 
discussion  given  here  draws  upon  both  of  these  works. 

In  describing  the  turbulent  or  nonturbulent  stream,  the 
velocity  and  pressure  are  separated  into  a  mean  value  component  and 
an  oscillating  component.  Instantaneous  velocity  and  pressure  values 
are  represented  as  the  sum  of  the  two  components.  The  mean  value 
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components  are  generally  defined  with  respect  to  time  at  a  fixed  point 
In  space.  Due  to  the  statistical  nature  of  turbulence,  they  may  also 
be  averaged  at  different  points  in  space  at  the  same  time.  The  total 
velocity  components  and  pressure  are  then  given  by  the  following: 


u  «=  U  + 

u' 

(1) 

v  ■  V  + 

v' 

(2) 

w  ■  W  + 

w' 

(3) 

p  ■  P  4 

P* 

(A) 

where  capitals  denote  the  mean  values  and  primes  the  fluctuating 
values.  Two-dimensional  mean  flow  with  two-dimensional  disturbances 
la  considered.  Schlichting  and  Schubauer  note  that  this  is  the  conserva 
tive  case  in  that  analysis  of  three-dimensional  disturbances  will  give 
higher  critical  Reynolds  Numbers.  Simplifications  are  introduced  by 
stating  that  V  *  0  everywhere  and,  in  the  boundary  layer,  the  mean 
flow  is  only  a  function  of  distance  from  the  surface,  y.  Thus,  we  have 


U  -  f(y) 

(5) 

V  *  0 

(6) 

w  *  0 

(7) 

i'  «  u'(x,y,t) 

(8) 

r'  *  v '  (x,y , t ) 

(9) 

»'  ■  P'(x,y,t). 

(10) 

It  is  assumeJ  that  both  the  total  motion,  given  by  Equations  1  through 
1,  and  the  mean  motion  are  solutions  of  the  Navier-Stokes  equations. 
Substituting  Equations  1  through  3  in  the  Navier-Stokes  equations  gives 
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the  following  two  equations  after  squares  and  products  of  disturbance 
velocities  and  of  their  differential  coefficients  are  neglected. 


3ul 

dt 


+  v 


,au 

by 


,2u')  •  Kf£ +  fir)  • 


an 


dv' 

at 


+  u 


av' 

ax 


vV  v' 


(12) 


The  continuity  equation  for  the  disturbances  is 


&L  + 

ax  ay 


o. 


(13) 


Since  the  mean  flow  satisfies  the  Navier-Stokes  equation,  two  more 
relations  are  gained  by  letting  u1  =  v'  «  p'  ■  0  in  Equations  11  and  12. 


0 


0 


(14) 


(15) 


Subtracting  from  Equations  11  and  12,  respectively,  Equations  13  and 
14,  gives  expressions  solely  in  terms  of  U  and  the  disturbances. 


at 


v72u' 


I  gP_l 

p  ax 


(16) 


av' 

at 


+  u 


av' 

ax 


vV  v' 


lb2l 
p  3y 


(17) 


Equations  16  and  17  arc  the  fundamental  hydrodynamic  equations  for 
snail  disturbances.  Those  equations  may  be  subtracted  to  eliminate 
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the  pressure  disturbance,  after  selective  differentiation,  yielding, 
together  with  the  continuity  equation,  two  equations  for  u'  and  v'. 

The  mathematical  form  of  the  disturbance  is  suggested  by  the 
nature  of  the  oscillations.  Observing  the  oscillations  at  a  fixed 
point  suggests  the  passage  of  a  traveling  wave  in  the  boundary  layer. 
Thus,  the  disturbance  is  periodic  and  may,  therefore,  be  represented 
by  a  Fourier  series  and  investigated  using  a  single  term  of  the  series. 
A  stream  function,  i{i,  which  satisfies  this  condition  is 

♦  ■  F(y)  exp[i(ax  -  0C )]  .  (18) 

The  amplitude  of  the  stream  function,  F,  is  represented  by  F(y). 

Since  the  mean  flow  depends  solely  on  y,  the  amplitude  is  also  assumed 
to  depend  on  y  only.  The  wavelength  is  given  by 

X  ■  2tt fa  ,  (19) 

and  the  frequency  and  amplification  are  given  by  p.  The  quantity  8  is 
complex 

e  -  er  ♦  ,  (20) 

where  p^  is  the  circular  frequency  or  2rrf,  and  p^  is  the  amplification 
coefficient.  By  substituting  Equation  20  in  Equation  18,  it  is  seen 
that  the  disturbances  grow  with  time  if  3^  is  greater  than  zero.  The 
wave  velocity,  c  ,  is  then 


or 


(21) 
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The  stream  function  gives  the  disturbance  velocities  as  follows: 


i  M. 
u  *  .  . 

By 

(22) 

i  B& 

v'  *  -  r1  . 
dx 

(23) 

u'  *  F*  exp[i(crx  -  Bt)]  , 

(24) 

v'  ■  -  iaF  exp[i(ofx  -  BO]  , 

(25) 

where  the  prime  on  F  denotes  differentiation  with  respect  to  y. 

It  is  immediate'  seen  that  the  continuity  equation  is 

satisfied.  After  p'  is  eliminated  from  Equations  16  and  17,  the 

disturbance  velocities  are  introduced  and  the  resulting  equation  is 

nondimensionalized  with  respect  to  the  maximum  free  strean  velocity, 

U  ,  and  the  boundary-layer  thickness,  6.  The  resulting  fourth  order 
m 

nonlinear  equation,  the  Orr-Sommerfeld  equation,  is  called  the 
stability  equation. 


(U-c) (F"  -  a2F)  -  U"F  -  (F""  -  2o2F"  +  o4F)  . 

6 


(26) 


Here 


B 

c*&*c  +  i  c ,  , 

a  r  i  ' 


Re 


U  6 
m 

6  “  v  * 
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The  terms  on  the  left  side  of  the  equation  represent  inertial  terms 
and  those  on  the  right  side  viscous  terms.  For  large  Reynolds  Number, 
the  viscous  terms  are  neglected  and  the  resulting  equation  is  known  as 
the  frictionless  stability  equation.  The  boundary  conditions  applicable 
to  Equation  26  demand  that  the  components  of  the  turbulent  velocity  be 
zero  at  the  wall  and  at  a  large  distance  from  it.  These  boundary  con¬ 
ditions  are  the  following: 

At  y  ■  0  and  y  *  ®; 

u '  ■  v 1  ■  F  «  F*  ■  0. 

If  the  frictionless  stability  equation  is  Investigated,  only  two  of 
the  above  boundary  conditions  may  be  satisfied. 

A  general  solution  of  Equation  26  may  be  formed  with  four 
particular  solutions.  The  form  of  the  particular  solutions  depends  on 
the  form  of  the  boundary-layer  velocity  profile.  Both  Tollmein  and 
Schllchting  obtained  solutions  for  profiles  closely  approximating  a 
Blaslus  dis  ribution.  The  solutions  reduce  to  an  eigen-value  problem 
in  terms  of  the  parameters  a ,  Re^,  and  c.  For  each  a  and  Re  value,  one 
complex  eigen-value,  c,  is  found.  The  locus  of  points  where  c^,  or 
essentially  equals  zero  is  the  curve  of  neutral  stability.  En¬ 
closed  by  this  curve  are  unstable  solutions.  Theoretical  wave  param¬ 
eters  for  neutral  oscillations  taken  directly  from  Schllchting  are 
given  below  in  Table  2.  Here  Re  is  based  on  the  displacement  thick- 
ness,  6  .  Using  the  tabulated  values,  Schlichting's  neutral  disturbance 
curve  is  plotted  in  terms  of  frequency,  wavelength,  and  velocity  in 
Figu  res  1  through  3. 
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TABLE  2.  NEUTRAL-STABILITY  CURVE 


Brai.ch  I _  _  Branch  II 


c 

r 

U 

m 

V 

* 

or  6 

0V  6 

iTT*  10 

m 

V 

* 

or6 

8  v 

fi*  10 

m 

0.200 

7200 

0.077 

2.14 

37,600 

0.149 

0.70 

0.250 

3010 

0.101 

8.39 

12,000 

0.188 

3.92 

0.300 

1530 

0.129 

25.3 

4,640 

0.223 

14.4 

0.325 

1150 

0.143 

40.4 

3,290 

0.238 

23.5 

0.350 

893 

0.159 

62.3 

2,070 

0.251 

42.4 

0.375 

736 

0.181 

92.2 

1,420 

0.264 

69.7 

0.400 

633 

0.205 

129.5 

1,020 

0.274 

107.4 

0.420 

606 

0.239 

165.7 

713 

0.273 

161.0 
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The  above  solution  for  neutral  oscillations  was  based  on  the 
assumption  that  U  c  f(y)  and  V  =  0  in  the  boundary  layer.  This  is 
equivalent  to  assuming  a  nonthickening  boundary  layer.  Since  the 
boundary  layer  grows  in  thickness  in  the  downstream  direction,  the 
validity  of  the  solution  may  be  questioned  on  physical  grounds.  This 
is  one  of  the  basic  criticisms  of  stability  theory.  However,  the 
dependence  of  U  on  y  is  large,  compared  to  the  dependence  on  x,  and  V 
is  small  compared  with  U.  The  experimental  work  of  Schubauer  and 
Skramstad  (1947)  qualitatively  confirmed  the  general  characteristics 
predicted  by  the  theory.  This  is  illustrated  in  Figure  1  where,  along 
with  Schlichting's  neutral-stability  curve,  are  plotted  experimental 
data  points  generated  by  ^chubauer  and  Skramstad.  The  qualitative 

* 

differences  are  partially  explained  by  the  fact  that  the  values  cf  6 
used  by  Schubauer  and  Skramstad  are  calculated  for  a  thickening  boundary 
layer  and  Increase  with  x.  Following  a  different  procedure,  developed 

g 

by  C.  C.  Lin  ,  Shen  repeated  the  calculations  for  the  neutral-stability 

curve  and  obtained  results  in  closer  agreement  with  experiment.  Shcn's 

curve  is  also  shown  in  Figure  1. 

For  the  present  investigation,  the  range  of  Re  *  is  suffi- 

6 

ciently  high  that  Schlichting's  results  are  quantitatively  satisfactory. 

★ 

Realizing  this  fact,  the  values  of  6  may  be  related  to  x  with  minimal 
compromise  of  theory  to  experiment.  For  convenience,  the  Reynolds 
Number  based  on  6  is  converted  to  Rc^  using  the  following  relation, 

6*  -  1.72  Qp)  7  ,  (27) 

m 


from  Reference  5. 
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Together  with 


U  x 

Re  =  ~ —  ,  and 
x  v 


Um  6 

Re  /.  *  -  this  gives 

6  v  ’ 


(28) 


Equation  28  was  used  to  replot  Figures  1  through  3  based  on  Re^  Instead 
of  Re^*.  The  resulting  curves  are  presented  in  Figures  4,  5,  and  6 
where  and  a  have  been  replaced  by  f  and  X,  respectively. 

Disturbances  that  enter  or  originate  within  the  bounds  of 
the  neutral-stability  curve  are  amplified.  This  amplification  process 
is  the  link  between  instability  and  transition.  The  amplification  of 
the  disturbance  may  be  expressed  as  the  ratio  of  the  amplitudes  at  two 
discrete  times.  From  Equation  18,  the  amplitude  ratio  at  two  consecu¬ 
tive  crests,  on  a  given  streamline,  is 


or 


A 

~  «=  exp(3it2)/exp(8it1)  , 


(29) 


expj  2  3 
C1 


dt 


) 


(30) 


where  is  the  amplitude  at  t^  and  is  the  amplitude  at  t^.  As 
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pointed  out  by  Gyorgyfalvy  ,  transition  can  be  predicted  empirically 
based  on  the  integrated  amplification  rate.  According  to  this 
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criterion,  transition  occurs  after  the  onset  of  instability  when 

/tr  dt  «=  9.  (31) 

For  disturbances  entering  at  Branch  I  of  the  neutral-stability  curve, 
Figure  1,  this  integrated  amplification  is  first  reached  at  Re^  «=  3.1 
x  10^.  This  is  in  good  agreement  with  experimental  values  for  transi¬ 
tion  on  a  flat  plate  with  zero  pressure  gradient  for  the  condition  of 
free  stream  turbulence  intensity  less  than  0.1,  Reference  5. 

The  theory  of  boundary-layer  oscillations  described  above 
forma  the  basis  for  this  experimental  investigation.  If  the  necessary 
condition  for  the  occurrence  of  tranaition  is  the  amplification  of 
boundary-layer  disturbances,  then  the  disruption  of  this  amplification 
process  may  be  expected  to  retard  transition.  The  possibility  of 
accomplishing  this  with  controlled-f requency  electrical  disturbances 
is  described  in  the  following  section. 
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ELECTRICAL  PERTURBATION’S 


The  sensitive  phenomenon  of  transition  may  be  influenced  by 

mechanical  disturbances  with  relatively  small  amounts  of  input  energy 

required.  This  is  illustrated  in  References  5  and  6  which  describe 

the  effects  of  loudspeakers  on  flow  transition.  In  addition  to  such 

10 

sonic  interference  investigations,  Deimen  and  Clark  have  experi¬ 
mentally  studied  the  effect  of  replacing  a  section  of  a  solid  flow 
boundary  with  a  vibrating  membrane.  All  of  the  above  efforts  have 
demonstrated  that  transition  may  be  promoted  by  mechanical  disturbances 
of  the  proper  frequency.  The  question  exists  as  to  what  other  forms 
of  disturbance  may  influence  transition  and  is  it  possible  to  forestall 
it. 

For  years  it  has  been  known  that  electric  fields  and  the 
presence  of  ions  effect  fluid  properties  and  performance.**  A  dramatic 
illustration  of  the  capability  of  corona  discharge  to  influence  boundary- 
layer  phenomena  is  given  in  Reference  12.  This  work  describes  the 
interaction  of  a  corona  point  discharge  with  a  free-convection  boundary 
layer  on  a  vertical  plate.  Significantly,  large  changes  in  heat- 
transfer  coefficient,  local  fourfold  increases  in  some  cases,  were  ob¬ 
tainable  with  corona  power  expenditures  of  less  than  1  watt.  A  recent 
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experimental  investigation  has  demonstrated  that  flow  attachment  may 
also  be  influenced  by  corona  wind.  Based  on  this  knowledge,  one  may 
conjecture  that  electrical  disturbances  may  be  capable  of  influencing 
transition.  For  this  study,  the  corona  discharge  was  chosen  as  the 
disturbance  oi  t cr turbot ion  mechanism. 


Corona  discharge  is  usually  viewed  as  undesirable.  When 
speaking  of  high-voltage  electrical  equipment  and  power  transmission, 
this  view  is  well  taken.  As  pointed  out  by  AIEE  engineers^,  corona 
power  loss  of  3  kw  per  mile  on  high-voltage  transmission  lines,  with 
accompanying  insulation  deterioration,  is  not  uncommon.  However,  many 
uses  of  corona  have  been  discovered  with  beneficial  results.  Several 
applications  are  air-pollution  control,  radiation  detection,  and  static 
elimination.  The  study  of  corona  discharge  has  also  served  science  in 
helping  to  explain  processes  of  electrical  breakdown  in  gases.  It  was 
the  positive-point  corona  which  led  physicists  to  the  discovery  of 
the  role  of  photoionization  in  gases.  A  description  of  the  corona 
discharge  and  its  application  to  the  experiment  follows. 

Corona,  as  defined  by  Loeb,  is  the  expression  used  to  describe 
the  class  of  luminous  phenomena  associated  with  the  current  jump  at  a 
highly  stressed  electrode  preceding  a  spark  breakdown  of  the  gap. 
Associated  with  corona  breakdown  is  an  electrode  of  small  radius  or 
high  curvature.  If  such  an  electrode  is  placed  near  one  of  low  curva¬ 
ture,  the  electric  field  induced  is  asymmetric  and  extremely  concen¬ 
trated  at  the  highly  curved  electrode.  As  the  potential  difference 
between  the  electrodes  is  raised,  ionization  of  the  air  in  the  immedi¬ 
ate  vicinity  of  the  highly  curved  electrode  occurs;  however,  no  spark 
crosses  the  gap  because  the  field  is  too  low  farther  from  the  curved 
electrode.  The  resulting  corona  discharge  will  be  of  positive  or 
negative  character  depending  upon  the  electrode  polarity. 
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If  the  highly  curved  electrode  or  wire  has  a  negative  po¬ 
tential,  the  positive  ions  formed  near  the  wire  will  acquire  relatively 
high  energy  from  the  field  and  essentially  bombard  the  cathode  wire. 
These  positive  ions,  whi le  forming  a  space  charge  near  the  cathode 
which  gives  the  negative  corona  its  fluctuating  characteristic,  will 
produce  the  necessary  electrons  to  sustain  the  discharge.  On  the  other 
hand,  the  electrons  that  are  formed  will  migrate  relatively  slowly  to 
the  anode  because  of  the  reduction  in  field  strength  away  from  the 
wire.  These  electrons  will  generally  not  have  enough  energy  to  cause 
further  ionization  and  may  be  expected  to  attach  themselves  to  neutral 
molecules  of  electro-negative  gases  forming  large,  even  slower  moving, 
ions.  A  wire  of  positive  potential  will  attract  the  ionized  electrons 
with  great  intensity  causing  electron  avalanches  toward  the  wire.  This 
liberates  a  positive  ion  space  charge  which  drifts  toward  the  cathode 
sustaining  the  discharge  although  reducing  the  field  at  the  wire. 

Positive  wire  corona  was  chosen  for  the  investigation  because 
of  the  more  uniform  nature  of  the  discharge.  It  was  feared  that  the 
nonperiodic  pulsating  nature  of  a  negative  wire  corona^  might  detri¬ 
mentally  effect  the  flow  due  to  noise  in  addition  to  adding  an  uncon¬ 
trolled  variable  to  the  experiment.  Qualitatively,  seme  relations 
between  the  corona  variables  arc  given  in  Table  3. 


16 


TABLE  3.  CORONA  VARIABLES 


Variable 

- .= 

Effect  of  Variable  Increase 
on  Corona  Current 

Gap  potential 

Increase 

Anode  curvature 

Increase 

Gap  spacing 

Decrease 

Density 

Decrease 

Once  liberated,  the  negatively  charged  bodies,  whether  In 
the  form  of  electrons  or  molecules,  will  experience  a  body  force  due 
to  the  presence  of  the  electric  field.  This  body  force  will,  in 
effect,  be  distributed  throughout  the  fluid  between  the  electrodes  by 
virtue  of  electron-molecule  and  ion-molecule  collisions.  In  this 
respect,  the  electrical  disturbance  is  physically  different  from  the 
mechanical  disturbances  mentioned  above. 

A  corona  discharge  with  a  controlled  frequency  component 


was  employed  with  the  aim  of  stabilizing  the  boundary  layer  as  described 
In  the  following  section.  Positive  ions  were  used  to  interact  with  the 
boundary-layer  flow. 
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EXPERIMENTAL  I  NAT  ST  I  GAT  I  ON 

DESIGN  OF  THE  EXPERIMENT 

In  attempting  to  interact  the  corona  discharge  with  the 
boundary-layer  instabilities,  an  essential  feature  of  the  experimental 
investigation  is  the  manner  in  which  the  electrical  disturbance  is 
applied  to  the  boundary-layer  flow.  Several  desirable  character¬ 
istics  are: 

(1)  The  disturbance  shall  incorporate  a 
pulsating  controlled-f requency  component. 

(2)  The  mechanism  for  applying  the  disturbance 
shall  not  trip  the  boundary  layer  by  virtue 
of  its  presence. 

(3)  The  disturbance  shall  cover  a  relatively 
large  area  of  the  flow  field. 

A  series  of  corona  wires,  suspended  outside  the  boundary  layer,  parallel 
to  the  flow  surface  and  normal  to  the  flow  direction,  appears  to  satisfy 
these  requirements. 

Several  variables  are  of  prime  importance.  These  are,  the 
disturbance  frequency  and  intensity  and  the  spacing  and  location  of  the 
wires  in  the  flow  direction.  Because  the  amplification  of  boundary- 
layer  oscillations  is  a  selective  process,  depending  on  frequency,  it 
would  seem  that  the  predominant  variable  of  the  disturbance  should  be 
the  disturbance  frequency.  However,  since  the  boundary-layer  oscilla¬ 
tions  are  in  the  form  of  waves,  cne  might  expect  that  any  disturbance 
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which  could  interact  with  the  oscillations  would  also  have  wavelike 
characteristics.  This  leads  to  the  idea  of  adjusting  the  corona  wire 
spacing  to  be  a  multiple  of  the  wavelength  which  is  possessed  by  a 
predominant  boundary-layer  oscillation.  The  provision  for  attaining 
this  spacing  is  described  in  a  later  section  along  with  the  electrical 
equipment  which  controls  the  corona  frequency  and  intensity. 

The  existence  of  predominant  oscillations  has  been  theorized 
by  Schlichting  and  demonstrated  by  Schubauer  and  Skramstad,  Reference  5, 
and  Deimen  and  Clark,  Reference  10.  Such  oscillations  obtain  because 
of  the  nature  of  the  amplification  process.  Referring  to  Figure  4, 
waves  of  many  frequencies  enter  the  amplification  zone  through  Branch  I 
and  are  then  amplified  at  a  rate  dependent  on  their  frequency.  They 
then  proceed  to  the  value  of  Re^  at  which  an  observation  is  made. 
Initial  oscillations  arc  supposedly  present  in  the  boundary  layer  at 
all  values  of  Re^  and,  therefore,  enter  the  diagram  at  all  points.  It 
is  obvious  from  Equation  30  that  oscillations  which  enter  at  Branch  I 
and  proceed  to  Branch  II  will  receive  the  highest  amplification.  When 
an  observation  is  made  at  a  specific  value  of  Re^,  there  exists  one 

value  of  $  v/U  2  which  will  be  more  highly  amplified  than  any  other. 

r  m 

This  is  the  predominant  oscillation.  The  question  that  remains  is 
where  on  the  diagram  is  the  locus  of  the  predominant,  or  most  ampli¬ 
fied,  oscillations.  This  can  be  determined  by  considering  the  initial 
energy  distribution  of  the  oscillations  and  the  damping  prior  to  reach¬ 
ing  Branch  I  of  the  neutral-stability  curved  The  desired  locus  of 
points  lies  in  the  amplification  zone  near  Branch  II.  Knowing  this, 


19 


the  corona  discharge  frequency  and  wire  spacing  may  be  gauged  to 
match  the  predominant  oscillations. 

Further  details  concerning  the  design  of  the  experiment  and 
apparatus  are  described  in  the  section  below. 

EXPERIMENTAL  APPARATUS 


The  experimental  apparatus  for  this  transition  study  was 
selected  or  designed  to  yield  the  following  test  conditions: 

Boundary  surface  -  flat  with  zero  pressure  gradient, 

electrical  conductor 

Flow  conditions  -  steady-state,  incompressible 

flow  of  air 

Transition  detector  -  reliable,  rapid  nonconductor 

Electrical  equipment  -  high  voltage,  low  current,  con¬ 
trolled  frequency  pulsation. 

Boundary  Surface 

The  boundary  surface  was  a  3-foot -wide  by  6-foot-long  flat 
plate  consisting  of  a  thin  sheet  of  aluminum  laminated  to  a  3/4-inch- 
thick  sheet  of  A-A  exterior  grade  plywood.  The  aluminum  was  3003  H14 
cold-rolled  sheet  1/32-inch  thick  and  36  inches  wide.  Surface  roughness 
values  for  such  cold-rolled  products  typically  range  from  10  to  20 
microinches,  rms.  Special  care  was  taken  to  keep  the  plate  surface 
scratch-frec  and  clean.  Waviness  was  imperceptible  even  when  viewed 
near  grazing  incidence;  however,  some  slight  sag  was  present  as  revealed 
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by  a  crack  of  light  showing  beneath  a  straight  edge  placed  on  the 
plate  normal  to  the  flow  direction.  A  straight  edge  placed  parallel 
to  the  flow  direction  revealed  no  cracks  of  light. 

The  leading  edge  of  the  plate  was  a  15-degrec  bevel. 

Initially,  the  bevel  nose  was  roundedj  however,  early  test  results 
led  to  the  subsequent  application  of  a  sharp  builtup  nose.  This  is 
described  in  the  testing  section.  Full-scale  drawings  of  both  leading 
edge  shapes  appear  in  Figure  7.  Three  rows  of  holes  for  static  pressure 
taps  were  placed  in  the  plate  surface  to  allow  determination  of  pressure 
gradients  on  the  plate.  A  layout  of  the  holes  and  a  section  view  of  a 
pressure  tap  is  shown  in  Figure  8.  A  curved  aluminum  flap  was  attached 
to  the  trailing  edge  of  the  plate  to  maintain  the  flow  blockage  neces- 
•ary  for  keeping  the  stagnation  point  on  the  active  side  of  the  plate. 
This  flap  may  be  seen  in  Figures  9  and  10. 

The  plate  was  mounted  horizontally  at  mid-height  in  the  flow 
channel  test  section  with  the  active  side  facing  downward.  Steel  angles 
were  used  to  secure  the  plate  in  the  test  section  with  a  minimum  of  flow 
interference,  Figure  11.  Slots  in  the  larger  angles  allowed  the  plate 
to  be  inclined  to  the  flow.  No  vibration  of  the  plate  was  visible 
during  testing. 


Flow  Channel 

The  flow  channel  used  for  this  experimental  investigation  is 
Battelle  Memorial  Institute's  subsonic  wind  tunnel,  the  general  layout 
of  which  is  shown  in  Figure  12.  The  tunnel  is  an  open-circuit,  or 
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"NPL"  type  which  draws  in  fresh  outside  air  continuously.  From  the 
settling  chamber  to  the  test  section,  the  area  reduction  is  approxi¬ 
mately  10:1.  A  125-hp  induction  motor  drives  the  fan  at  700  rpm. 

The  fan  (Figure  13)  is  a  Hartzell  No.  MP-96-6  with  a  diameter  of  8 
feet  and  adjustable  pitch  blades.  All  testing  was  done  with  a  blade 
pitch  of  17  degrees.  The  velocity  profile  across  the  tunnel  test 
section  is  flat.  This  can  be  seen  from  a  typical  plot  of  dynamic 
pressure  ratio  versus  distance,  Figure  14.  In  the  test  section,  the 
wall  boundary  layer  is  turbulent  and  several  inches  thick.  Sec 
Figure  15. 

The  level  of  turbulence  in  the  test  section  is  low,  as 
governed  by  three  24-mesh  screens,  of  0.0075-inch  wire,  installed  at 
the  entrance  to  the  settling  chamber.  The  projected  open  area  of  the 
screens  is  67.4  percent.  Figures  16  and  17  show  the  tunnel  mouth  and 
turbulence  screens.  The  method  employed  to  measure  the  turbulence 
made  use  of  a  5-inch-diameter  "pressure  sphere".  Pressure  sphere 
measurements  yield  the  turbulence  factor  which  is  defined  as  the  ratio 
of  an  effective  Reynolds  Number  to  the  actual  Reynolds  Number.  For 
turbulence  factors  less  than  2,  the  percent  turbulence  is  related  by 
the  following  expression^. 


Percent  turbulence 


Percent  turbulence 


1.25  (turbulence  factor-1). 

2  2  2  1/2 
[  (u ' )  +  (v')  +  (w')  ]  x  100/U  . 

in 


(32) 

(33) 
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At  a  flow  speed  of  140  fps,  pressure  sphere  determinations 
yielded  a  turbulence  factor  of  1.083.  This  relates  to  a  turbulence 
intensity  of  0.10  percent.  In  general,  the  turbulence  level  increases 
with  tunnel  test  velocity. 

Transition  Detector 

Many  methods  may  be  used  to  detect  the  locations  of  the 
transition  region.  Some  of  the  most  common  include  the  following: 

(1)  Measuring  or  plotting  the  velocity  profile  in 
the  boundary  layer  and  determining  whether 
the  flow  is  laminar  or  turbulent  by  the 
profile's  shape. 

(2)  Reading  the  dynamic  pressure  at  a  small 
height  above  the  surface  and  noting  the 
distance  from  the  leading  edge  where  it  is 
a  minimum. 

(3)  Reading  the  dynamic  pressure  as  in  (2),  or 
velocity,  with  a  hot-wire  anemometer,  and 
noting  unsteady  or  oscillating  values. 

(4)  Emitting  smoke  from  the  surface  and  noting 
the  point  where  the  smoke  is  dispersed. 


*  Personal  communication  with  V.  0.  Hochne,  Battelle  Memorial 
Institute,  Columbus,  Ohio. 
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(5)  Coating  the  flow  surface  with  evaporating 
chemicals.  Where  the  flow  is  turbulent, 
the  evaporation  will  proceed  most  rapidly. 

(6)  Listening  to  a  total  head  probe  with  a 
stethoscope  as  the  probe  is  moved  along 
the  flow  surface.  A  laminable  flow  makes 
a  hissing  sound  while  a  roar  is  heard  when 
the  flow  is  turbulent. 

A  near-surface  impact  pressure  probe  was  used  to  detect 
transition,  making  use  of  Methods  (2)  and  (3)  above.  This  method  was 
chosen  because  it  is  simple,  rapid,  visual,  and  amenable  to  data 
recording. 

If  an  impact  pressure  probe  is  traversed  along  the  flow 
surface,  in  the  flow  direction,  the  reading,  which  is  directly  pro¬ 
portional  to  the  square  of  the  local  velocity,  is  determined  by  the 
thickness  of  the  boundary  layer,  the  height  of  tne  probe  above  the 
surface,  and  the  shape  of  the  local  velocity  profile.  As  the  probe 
moves  away  from  the  leading  edge,  at  a  fixed  distance  from  the  surface, 
the  laminar  boundary  layer  thickens,  effectively  moving  the  probe  into 
abions  of  lower  velocity.  However,  once  the  turbulent  tone  is 
readied,  the  velocity  profile  changes,  causing  an  abrupt  increase  in 
local  velocity  and  dynamic  pressure.  As  the  probe  continues  moving 
rearward  in  the  turbulent  region,  the  dynamic  pressure  again  decreases 
due  to  the  thickening  of  the  turbulent  boundary  layer.  This  sequence 
of  events  is  illustrated  in  Figure  18.  No  definite,  fixed  transition 
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point  exists.  The  start  of  the  transition  region  is  taken  as  the 
minimum  point  on  the  curve,  Point  (a).  Additionally,  this  point  is 
marked  experimentally  by  oscillations  in  velocity  readings  because  of 
the  unsteady  nature  of  the  transition  process. 

A  motorized  probe  was  designed  to  traverse  the  underside  of 
the  plate  in  the  longitudinal  direction.  The  probe  boom,  a  1-1/4- 
inch-square  aluminum  extrusion,  was  suspended  and  pivoted  on  a  crossbar 
downstream  of  the  plate  and  test  section.  The  crossbar  penetrated  the 
wind  tunnel  walls  and  was  isolated  from  tunnel  structure  vibrations. 

A  sprockct-and-chain  drive  was  employe^  capable  of  moving  the  boom  at 
a  rate  of  10  inches  per  minute.  The  probe  mechanism  may  be  seen  in 
Figures  9  and  10.  At  the  forward  end  of  the  boom,  an  offset  steel 
tube  accepted  the  glass  probe  used  for  the  measurements.  The  jlass 
portion  of  the  probe  was  connected  directly  to  a  pressure  line  which 
ran  through  the  boom  and  lead  out  of  the  tunnel. 

Figure  19  shows  a  typical  glass  probe  insert.  The  glass 
probe  tip  could  be  sized  to  virtually  any  length  and  diameter  by 
suitable  torch  heating  and  drawing.  Tip  diameters  ranging  from  0.010 
to  0.050  inch  were  used  in  the  experiment.  The  tip  ends  were  ground 
flat  using  a  small  croakus  cloth  wheel  in  an  electric  drill.  In 
Figure  20,  the  probe  offset  and  glass  insert  are  shown  in  relation  to 
the  plate  surface.  The  offset  was  adjusted  to  maintain  the  glass  tip 
tangent  to  the  plate  at  all  times. 
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Electrical  Equipment 

A  controlled  frequency  pulsating  potential  was  applied  to  a 
series  of  corona  wires  suspended  over  the  flat  plate.  The  wires  were 
suspended  and  tightened  by  means  of  a  double  row  of  screws  which  were 
installed  in  plexiglass  strips,  Figures  20  and  21,  secured  to  the 
underside  of  the  plate-suppor t  angles.  Within  the  test  section,  all 
exposed  conductors  were  grounded  with  the  exception  of  the  corona 
wires  and  rows  of  screws.  Several  types  of  corona  wires  were  tried 
during  the  course  of  the  experiment.  Ultimately,  the  wire  most  resis- 
tant  to  breaking  proved  to  be  0.008-inch  steel  model  airplane  control 
wire.  Music  wire  cf  0.004-inch  and  0.011-inch  diameter  proved  to  be 
unsatisfactory. 

A  pulsating  potential  was  applied  to  the  wires  by  direct 
superposition  of  an  alternating  voltage  signal  on  a  steady  d-c  voltage. 
This  input  form  is  illustrated  in  Figure  22.  The  superposition  was 
employed  because  of  the  difficulty  and  expense  involved  in  obtaining 
or  designing  equipment  to  yield  a  pulsating  signal  having  a  potential 
level  as  high  as  the  d-c  level.  A  circuit  diagram  appears  in  Figure  23. 
The  components  of  the  electrical  apparatus  may  be  seen  in  Figure  24  and 
are  specified  below. 

Oscillator  -  Hewlett-Packard  wide-range  oscillator,  Model  200  CD. 

This  generates  the  variable  frequency  component  of 
the  corona  potential. 
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Amplifier  -  Scott,  100-watt,  laboratory  power  amplifier, 

Tvpc  265A. 

The  amplifier  boosts  the  oscillator  output  for 
input  to  the  stepup  transformer. 


Transformer  -  Merit,  25-watt,  audio  frequency  transformer, 

Model  A-3129. 

This  transformer  was  run  "backwards"  to  gain  the 
maximum  voltage  stepup  possible. 


High  Voltage  D-C 
Power  Supply 


Custom  designed  for  Battelle  Memorial  Institute. 
Maximum  output  -  14,000  volts,  ripple-free. 


Voltage  Divider  -  Voltage  split  of  50,000:1,  total  resistance  of  50 

million  ohms. 

Used  in  conjunction  with  oscillograph  for  voltage 
measurement. 

Scope  -  Dumont  oscillograph,  Model  304A. 

Used  to  measure  and  visualize  imposed  voltage. 

Ammeter  -  Simpson,  Model  374  microammeter. 

Used  to  measure  average  corona  current. 


The  performance  of  the  electrical  system  varied  with  the 
size  and  number  of  corona  wires  used.  Normally,  four  0.008-inch  wires 
were  used,  each  36  inches  long.  With  such  a  load,  the  maximum  d-c 
output  attainable  was  12,000  volts.  The  maximum  pcak-to-peak  fluctu¬ 
ating  voltage  attainable  was  a  function  of  frequency  and  varied  from 
almost  nothing  at  1  cps  to  4000  volts  at  6000  cps. 
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EXPERIMENTAL  PROCEDURE 

Before  proceeding  with  transition  tests,  the  static  pressure 
distribution  over  the  surface  of  the  plate  was  checked,  in  the  presence 
of  the  traversing  pressure  probe,  at  several  flow  speeds.  The  surface 
pressure  taps  (see  Figure  8)  were  connected  to  the  manometer  board 
shown  in  Figure  25.  Each  U-tube  contained  Merriam  Red  Oil,  specific 
gravity  0.827,  and  was  open  to  the  atmosphere  on  one  side.  The  column 
height  could  be  read  conveniently  to  0.05  inch  which  is  the  limit  of 
accuracy  for  simple  manometers.  Both  legs  of  the  U-tube  were  read 
simultaneously.  Generally,  the  static  pressure  readings  were  within 
0.1  inch  of  the  average  of  the  readings,  the  maximum  deviation  being 
less  than  two  percent  of  the  total  dynamic  head.  Figure  26  presents 
one  set  of  typical  test  data  for  the  static  pressure  distribution  over 
the  plate.  The  static  pressure  readings  were  recorded  at  the  beginning 
of  each  test  and  remained  constant  throughout  the  test  for  a  given  speed 
setting. 

For  the  transition  tests,  the  experimental  procedure  is  out¬ 
lined  below. 

(1)  The  desired  tunnel  speed  was  set  approximately. 

When  equilibrium  temperature  and  flow  conditions 
were  reached,  the  speed  was  set  more  accurately. 

(2)  Transition  was  probed  for  using  the  traversing 
pressure  probe. 

(3)  The  desired  corona  wires  were  installed  followed 
by  a  reprobc  for  transition  detection. 
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(4)  The  corona  discharge  was  activated  followed 
by  a  third  probe  for  transition  detection. 

Details  relating  to  Parts  (1)  through  (4)  of  the  experi¬ 
mental  procedure  are  described  in  the  section  entitled,  RESULTS  OF 
EXPERIMENTAL  INVESTIGATION.  Some  of  the  experimental  calculations 
which  were  required  to  perform  the  tests  are  presented  in  the  follow 
ing  paragraphs. 


EXPERIMENTAL  CALCULATIONS 


To  perform  the  transition  tests  required  determination  of 
the  following  quantities: 

(1)  Air  density. 

(2)  Flow  speed. 

(3)  Local  Reynolds  Number. 

(4)  Frequency  and  wavelength  of  predominant 
boundary-layer  oscillations. 

Air  Density 

It  is  necessary  to  know  the  air  density  to  calculate  the 
flow  speed  and  the  local  Reynolds  Number.  The  tests  were  conducted 
at  flow  speeds  of  approximately  150  to  200  fps  at  temperatures  close 
to  30  degrees  F.  Corresponding  Mach  Numbers  are  less  than  0.2, 
classifying  the  flow  as  subsonic  and  incompressible.  The  air  density 
is  given  by  the  equation  of  state  for  a  perfect  gas. 
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P 


(34) 


For  convenience,  the  pressure  is  left  in  barometric  units. 
The  absolute  pressure,  p  ,  is  determined  by  subtracting  the  local 
static  pressure  near  the  plate  from  the  barometric  pressure. 


.lb  .  p  (in.  hg)  x  70.8(lb/ft^/in.  hg)  p 

P \T~3/  =  53. 3 (ft -lb, /lb  - °R)  x  T(°R)_  =  1,33  T  * 

1 1  i  m 


(35) 


The  air  temperature  within  the  test  section  is  displayed  at 
the  tunnel  control  console. 


Flow  Speed 

The  test  section  flow  speed  is  calculated  by  applying  an 
experimentally  determined  correction  factor  to  Bernoulli's  equation. 
Bernoulli's  equation  is 

p  u  2 

p  +  — jj —  •=  constant  .  (36) 


gives 


By  assuming  zero  velocity  in  the  settling  chamber,  this 


h 


q 


(37) 


where  h  is  the  difference  in  static  pressure  across  the  inlet  con- 

q 

traction.  Both  h  and  the  static  pressure  near  the  entrance  to  the 

q 

test  section,  h  ,  are  displayed  at  the  tunnel  control  consol  manom- 
s 

etors,  Figure  27.  The  flow  speed  correction  factor  accounts  for 
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boundary- lay er  growth,  frictional  effects  and  settling  chamber  veloc¬ 
ity.  Using  total  head  probes,  the  correction  factor  for  the  Battello 
wind  tunnel  has  been  found  to  be  1.01.  Thus,  Equation  37  is  inverted 
and  corrected  to  read, 


U 

m 


1.01 


1/2 


<38) 


In  terms  of  the  direct  reading,  h  ,  this  becomes 


or 


U 

m 


(39) 


The  tunnel  head,  h^,  was  varied  by  adjusting  the  position  of  the  air- 
intake  ports  shown  in  Figure  28.  These  openings  are  approximately  8 
feet  upstream  of  the  fan. 

Local  Reynolds  Number 

The  local  Reynolds  Number  is  defined  in  terms  of  the  distance 
from  the  leading  edge  of  the  plate. 


Re 


U  x 
m 


(40) 
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The  only  unprescribed  quantity  in  Equation  40  is  the  viscosity  p, .  The 

viscosity  was  considered  to  be  a  function  of  temperature  only.  A 

19 

linear  variation,  prescribed  by  the  following  points,  used 


Air  Viscosity,  lb^/ft-sec 
1.110  x  10'5 
1.285  x  10‘5 


Temperature,  °F 

0 

100 


Frequency  and  Wavelength  of  Boundary-Layer  Oscillation 

The  amplification  characteristics  of  boundary-layer  oscilla¬ 
tions  are  presented  in  Figures  4  and  5  in  terms  of  Re^  and  the  fre- 

2  * 

qucncy  and  wavelength  parameters  fv/U^  and  X/6  ,  respectively.  At  a 
given  Reynolds  Number,  it  is  desired  to  know  the  specific  wave  char¬ 
acteristics  for  selected  points  on  the  diagram.  The  frequency  is 

2  2 

obtained  directly  by  multiplying  the  values  fv/U  by  U  /v. 

m  m 

Equation  35  is  used  to  obtain  the  kinematic  viscosity.  The  wavelength 

•k 

is  determined  from  the  dimensions  parameter  \/5  by  using  Equation  27 
★ 

to  find  6  . 


32 


RESULTS  OF  EXPERIMENTAL  INVESTIGATION 

The  experimental  investigation  may  best  be  described  by 
three  categories  of  chronological  effort: 

(1)  Tests  and  modifications  to  obtain  the  neces¬ 
sary  laminar  boundary  layer  on  the  plate. 

(2)  Tests  and  modifications  to  enable  the  de¬ 
tection  of  transition. 

(3)  Tests  and  modifications  to  determine  the 
effects  of  corona  discharge  on  transition. 

An  account  of  these  tests  and  results  constitutes  the  remainder  of 
this  section. 

OBTAINING  A  LAMINAR  BOUNDARY  LAYER 

A  flow  region  with  an  associated,  stable  laminar  boundary 
layer  is  mandatory  for  investigating  transition  phenomena.  Several 
tests  were  run  before  it  was  determined  that  such  a  flow  region  did 
not  exist  with  the  initial  leading  edge  geometry  shown  in  Figure  7a. 
The  traversing  pressure  probe  was  used  with  glass  tip  diameters  of 
0.01  inch  and  0.04  inch  at  several  flow  speeds.  Traverses  were  con¬ 
ducted  along  the  plate  centerline  moving  from  the  rear  of  the  plate 
toward  the  leading  edge.  The  pressure-distance  pattern  shown  in 
Figure  18  was  sought.  The  pressure  dip,  or  pressure  fluctuation, 
signaling  the  onset  of  transition  was  expected  to  occur  between  1  and 
2  feet  from  the  leading  edge.  Instead,  pressure  fluctuations  were 
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observed  several  inches  from  the  leading  edge.  Further  back  from  the 
leading  edge,  no  transition  could  be  detected.  Data  illustrating 
these  fluctuations  are  presented  in  Figure  29  for  three  different 
tests.  The  erratic  pressure  distributions  near  the  leading  edge  indi¬ 
cated  the  presence  of  premature  transition,  or  possibly  separation. 

At  this  point  it  was  decided  to  modify  the  leading  edge  be¬ 
cause  it  was  felt  that  the  rounded  shape  was  causing  the  irregular 
flow.  A  sharp  edge  was  modeled  onto  the  existing  nose  section  using 
a  polyester  base  automobile  body  putty,  Figure  7b.  The  final  surface 
finish  was  obtained  by  alternately  buffing  with  triple-zero  wet  paper 
and  spraying  with  silver  paint.  The  surface  change  from  the  painted 
edge  to  the  aluminum  was  visually  apparent  but  barely  perceptible  to 
the  fingertips.  With  a  new  glass  probe  installed,  0.042-inch  diameter, 
several  tests  were  performed  to  again  seek  the  transition  impact 
pressure  pattern.  In  Tests  5  and  6  ,  the  plate  was  tilted  to  slight 
positive  and  negative  angles  of  attack.  Tilting  the  plate  caused  a 
perceptible  charge  in  the  static  pressure  distribution  In  the  flow 
direction,  showing  a  negative  gradient  for  +  1/2-degree  angle  of 
attack  and  a  positive  gradient  for  -  1/2-degree  angle  of  attack.  How¬ 
ever,  in  all  cases,  the  irregular  impact  pressure  pattern  near  the 
nose  was  still  present.  Figure  30  displays  the  pressure  profiles  for 
the  three  tests  with  the  sharp  leading  edge. 

It  is  believed  that  the  severe  dips  in  impact  pressure  shown 
in  Figure  30  are  caused  by  boundary-layer  separation  at  the  sharp 
leading  edge.  If  the  surface  probe  was  within  a  separation  bubble,  the 
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local  velocity  might  well  be  directed  opposite  to  the  flow  direction 
due  to  eddying.  This  would  explain  the  severe  pressure  dips.  Such 
separation  would  be  caused  if  the  flow  stagnation  point  rested  on  the 
inclined  edge  of  the  nose  rather  than  at  the  top.  The  flow  would  then 
have  to  turn  abruptly  around  the  sharp  leading  edge  in  order  to  pass 
to  the  active  side  of  the  plate. 

Another  modification  of  the  plate  was  then  undertaken  to 
remedy  this  condition.  A  curved  aluminum  flap  was  secured  to  the 
trailing  edge  of  the  plate  as  shown  in  Figures  9  and  10.  This  proved 
to  be  the  solution,  for  the  irregular  pressure  pattern  near  the  leading 
edge  was  no  longer  present.  Test  data  confirming  this  fact  are  also 
presented  in  Figure  30.  It  in  believed  that  the  flap  provided  suffi¬ 
cient  blockage  between  the  active  surface  of  the  plate  and  the  test 
section  floor  to  force  the  impinging  streamlines  to  curve  upward  when 
approaching  the  plate,  thus  forcing  the  stagnation  point  to  the  active 
side  of  the  plate.  The  presence  of  the  flap  did  not  distort  the 
static  pressure  distribution.  This  concludes  the  description  of  the 
test  efforts  concerned  with  obtaining  a  laminar  boundary  layer. 

DETECTING  TRANSITION 

Although  Test  7  (Figure  30)  verified  the  existence  of  a 
smooth  boundary-layer  development,  it  did  not  indicate  the  presence 
of  transition.  Another  test  was  conducted  with  the  same  0.042-in. - 
diameter  probe  at  a  flow  speed  of  165  fps.  No  transition  was 
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detected.  The  glass  probe  was  replaced  with  or.e  having  an  external 
tip  di  ameter  of  0.022  inch  with  a  0.011-inch  hole.  Using  this  probe 
impact  pressure,  fluctuations  were  observed  in  the  next  test,  Test  8, 
at  1.65  feet  from  the  leading  edge.  Figure  31  shows  the  pressure  pro¬ 
file  with  the  extreme  value  of  pressure  fluctuation  plotted  at  1.65 
feet.  The  pressure  fluctuations  were  persistent  with  a  peaking  fre¬ 
quency  of  approximately  once  per  second.  The  distance  from  the  leading 
edge  where  these  fluctuations  were  encountered  was  taken  as  the  start¬ 
ing  point  for  transition.  All  subsequent  tests  were  conducted  with  a 
probe  tip  diameter  equal  to  0.022  inch. 

Following  Test  8,  in  which  transition  was  detected,  the  tra¬ 
versing  impact  pressure  probe  was  ••balanced"  against  a  fixed  impact 
pressure  probe.  The  fixed  probe  was  mounted  as  shown  in  Figure  32. 

This  was  done  in  order  to  partially  eliminate  the  effects  of  tunnel 
gusts  on  windy  days.  Test  9  was  run  at  164  fps.  Transition  was  de¬ 
tected  at  1.5  feet  from  the  leading  edge.  The  Impact  pressure  profile 
is  also  shown  in  Figure  31.  Associated  with  Tests  8  and  9  are  the 
following  transition  Reynolds  Numbers: 

Test  8,  Re  .  *  1.66  x  10^. 

’  xtr 

Test  9,  Re  .  *  1.54  x  10^. 

According  to  hot-wire  studies  performed  for  a  flat  plate  with  zero 
pressure  gradient  (Reference  5),  the  free  stream  turbulence  intensity 
associated  with  the  above  transitions  is  about  0.3  percent.  Although 
quite  low,  this  val  "  is  three  times  the  0.1  percent  value  quoted 
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earlier  for  a  140  fps  flow  speed.  Possibly  the  sensitivity  of  the 
pressure  sphere  technique  for  detecting  extremely  low  turbulence 
levels  is  not  sufficient. 

At  this  time,  it  was  desired  to  confirm  the  presence  of 
transition  before  undertaking  the  corona  tests.  A  visualization 
technique  was  tried  in  which  a  subliming  solid  is  used  to  coat  the 
surface  of  the  plate.  In  the  turbulent  region,  the  solid  will  sublime 
first,  leaving  a  visual  pattern  of  transition.  Solutions  of  alcohol, 
saturated  with  both  naphthalene  and  biphenyl, were  used.  After  appli¬ 
cation  of  the  solution,  the  alcohol  evaporated  leaving  t tie  surface 
coated  with  fine  white  crystals.  The  sublimation  temperature  of  both 
chemicals  is  close  to  70  degrees  F.  A  test  was  ..onducted  with  each 
chemical  at  an  ambient  temperature  of  approximately  55  degrees  F.  In 
both  cases,  the  crystals  failed  to  sublime  preferentially.  Instead, 
much  of  the  coated  strip  remained,  even  at  the  conclusion  of  the  test. 
The  air  temperature  was  obviously  too  low  to  allow  valid  use  of  this 
technique. 

An  analytical  correlation  was  used  to  verify  the  presence  of 
transition.  The  basis  of  the  correlation  is  as  follows:  At  a  given 
distance  from  the  leading  edge  and  from  the  plate  surface,  the  ratio 
of  local  velocity  to  free  stream  velocity  may  be  computed  from  the 
impact  pressure  probe  reading.  This  velocity  may  be  compared  to  those 
expected  for  laminar  and  turbulent  velocity  profiles  at  the  particular 
point  of  measurement. 
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The  local  velocity  ratio,  U/U^,  is  equal  to  the  square  root 
of  the  local  dynamic  pressure  ratio. 


U 

m 


(Full  impact  pressure  -  local  impact  pressure^ 
Full  impact  pressure  / 


(41) 


The  thickness  of  local  laminar  and  turbulent  boundary  layers  is  given 
in  Reference  7  for  growth  starting  at  the  leading  edge. 


6 


6 


T 


0.37  x 


(42) 


(43) 


For  the  laminar  flow,  the  local  velocity  at  a  given  distance  from  the 
surface  is  computed  using  the  similarity  parameter,  T). 


11  *  5  ‘i'eff.ctlve'V- 


(44) 


The  parameter,  1),  is  related  to  the  laminar  velocity  ratio,  (U/U  )  , 

m  1> 

as  tabulated  below. 


TABLE  4.  VELOCITY  RATIOS  FOR  A  LAMINAR  BOUNDARY 
LAYER  ALONG  A  FLAT  PLATE 


Ti 

<U/Vi 

0.0 

0.0 

1.0 

0.33 

2.0 

0.63 

3.0 

0.85 

4.0 

0.96 

5.0 

1.00 
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Using  Equations  41  through  44,  the  results  for  Test  9  are 
cited.  In  measuring  impact  pressure  near  a  wall,  the  effective  center 
of  the  probe  is  displaced  away  from  the  surface  because  of  the  local 
transverse  velocity  gradient.  Young  and  Mass  have  found  experimentally, 
for  a  flat-edged  probe that 


§*  -  0.131  4  0.0083  ~  ,  (45) 

o  o 

where  Ay  is  the  upward  displacement,  and  dQ  and  d^  are  the  outer  and 
inner  diameters  of  the  probe.  For  Test  15, 

d  -  0.022  in.  «  0.00183  ft. 

o 

dt  «  0.011  in.  -  0.000915  ft. 


This  gives  Ay  •  0.000158  ft  and 


d 

y  ce  ■  "T  +  Ay  *  0.001073  ft. 

■'effective  2 

Calculations  were  performed  at  three  distances  from  the  leading  edge. 
The  results  are  noted  in  Table  5. 


TABLE  5.  TRANSITION  CONFIRMATION*,  164  FPS  FEE.  PER  SECOND 


Parameter 

Distance,  Feet 

1/2 

1-1/2 

4 

V  fcet 

0.00345 

0.0060 

0.00978 

6t,  feet 

0.0133 

0.032 

0.0704 

T) 

1.70 

1.00 

0.60 

<U/Ua,>T 

0.706 

0.623 

0.57 

<u/Vl 

0.55 

0.33 

0.20 

factual 

0.57 

0.35 

0.49 

*  Test  9. 
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It  is  seen  tha'  up  to  the  transition  point,  the  actual 
velocity  rati  ;  matches  that  for  a  laminar  boundary  layer  quite  closely. 
The  boundary  layer  up  to  1.5  feet  is  clearly  not  turbulent  according 
to  the  calculations.  At  4  feet  from  the  leading  edge,  the  boundary 
layer  is  far  from  laminar  (0.49  versus  0.2)  showing  that  transition 
has  occurred.  On  the  basis  of  this  analytical  correlation,  the  ex¬ 
perimental  technique  for  detecting  transition  was  confirmed  and  tests 
involving  the  corona  apparatus  were  initiated. 


THE  EFFECT  OF  ELECTRICAL  PERTURBATIONS  ON  TRANSITION 

In  order  to  conduct  the  electrical  disturbance  tests,  or 
corona  tests,  it  was  necessary  to  select  a  flow  speed  and  an  associ¬ 
ated  disturbance  with  which  to  excite  the  boundary  layer.  At  this 
time,  it  was  decided  to  perform  all  subsequent  tests  at  a  flow  speed 
of  175  fps.  A  test  was  run  at  175  fps  and  transition  was  detected  at 
1.33  feet.  With  this  determined,  a  transition  Reynolds  Number  wa3 
estimated  as  follows: 


Re 


Barometric  Pressure,  29.5-in.  hg 
Temperature,  30*F 

Flow  Speed,  175  fps 

Transition  Distance,  1.33  feet 

0.078  x  175  x  1.33 


pU  x 
m 


xtr 


1.163  x  10’ 


1.56  x  10°. 


AO 


This  corresponds  to  Rc  *  *  2140. 

6  tr 


The  electrical  disturbance  was  chosen  to  match  the  fre¬ 
quency  and  wavelength  of  the  predominant  boundary-layer  oscillation 
at  this  Reynolds  Number  and  would  be  applied  to  the  flow  in  the  pre¬ 
transition  region.  As  discussed  above,  the  precise  location  of  the 
predominant  oscillation  is  not  known.  However,  it  is  known  that  the 
locus  of  predominant  oscillations  lies  close  to  Branch  II  of  the 

neutral  stability  curve  within  the  amplified  region.  A  point  was 

2r 
X 


ft  2tt  ★ 

chosen  for  a  value  of  or  6  *  T“  6  *  0.24.  Using  Equation  27, 


1/2 


.*  ,  „  /I -33  1.163  x  10  . 

6  ’  ll72  \  Vi  «  0.078  )  '  '’•°0185  £cct' 


This  gives  X  *  0.048  feet  and  X/6  *  26.6.  The  point  is  shown  on 

Figure  5.  To  determine  the  frequency  associated  with  the  wavelength 
chosen,  Figure  33  was  constructed  from  information  presented  by 

2 

Schubauer,  Reference  3.  It  is  seen  from  the  figure  that  g  v/U  * 

r  m 

37.4  x  10  ^at  Re  *  =  2140.  Thus,  the  excitation  frequency  is 

6 

0r/2n,  or 


37.4  x  10'*  U  2 


f  - 


2n  v 


m  37.4  x  10"6  x  1752  x  0.078 
—  ■  — —  1  ^ -  ■  1215  cps. 


2rr  x  1.163  x  10 


To  completely  characterize  the  oscillation,  the  wave  speed  is  computed. 


cr  -  f  X  -  1215  (0.048)  *  58.4  fps, 


c  /U  *=  0.334. 
r  m 


and 
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The  selected  disturbance  is  shown  as  a  point  on  Figures  4, 

5,  and  6  at  Re  =  Re  ;  the  maximum  total  amplification  associated 
with  a  wave  of  the  above  frequency  is  approximately  1000. 

In  order  to  isolate  electrical  effects,  it  was  first  neces¬ 
sary  to  conduct  several  "control"  tests.  The  screw-studded  plexiglass 
strips  were  installed  beneath  the  plate  and  a  test  was  conducted  at 
17S  fps.  The  location  of  transition  was  unaltered,  1.33  feet.  This 
is  shown  in  Figure  34.  Next  the  effects  of  wires  suspended  above  the 
boundary  layer  were  determined.  Two  possible  mechanisms  which  could 
trigger  a  premature  transition  were  present,  boundary-layer  entrain¬ 
ment  by  the  wake  behind  the  wires,  and  physical  vibration  of  the  wires 
Twenty-two  wires  were  installed,  at  a  longitudinal  spacing  of  0.048 
feet,  covering  the  region  from  0.917  to  1.940  feet  from  the  leading 
edge.  The  wires  were  0.004  inch  in  diameter.  Each  wire  was  spaced 
9/16  of  an  inch  from  the  surface  of  the  plate  to  allow  clearance  for 
the  1/2-inch-diameter  portion  of  the  pressure  probe.  This  9/16-inch 
distance  was  held  for  all  tests  involving  corona  wires.  Several  tests 
were  conducted  at  175  fps  with  the  wifes  installed.  Transition  was 
located  consistently  at  0.834  feet  from  the  leading  edge.  This  was 
immediately  in  front  of  the  first  wire.  It  is  believed  that  some 
distance  downstream  from  the  normal  transition  point,  the  rapidly 
thickening  boundary  layer  and  the  wake  left  by  the  wires  intermixed, 
entraining  the  boundary  layer  to  the  plane  of  the  wires.  This  entrain 
ment  propagated  upstream  until  the  first  wire  was  reached.  To  demon¬ 
strate  that  the  premature  transition  was,  in  fact,  caused  by  the  down¬ 
stream  wires,  one  0.008-inch  wire  was  suspended  at  a  distance  of 
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1.3  feet  from  the  leading  edge.  Transition  occurred  downstream  of  the 
wire  at  1.417  feet  from  the  leading  edge.  The  data  from  the  above 
tests  are  also  shown  on  Figure  34.  Additional  evidence  is  presented 
with  the  descriptions  of  the  corona  tests  substantiating  the  fact  that 
a  series  of  wires  may  be  suspended,  above  the  boundary  layer,  in  the 
pretransition  region  without  tripping  transition. 

The  corona  tests  were  initiated  by  suspending  three  wires, 
0.008  inch  in  diameter,  at  distances  of  1.209}  1.258,  and  1.307  feet 
from  the  leading  edge.  A  test  was  conducted  at  zero  current  for  which 
normal  transition  was  expected  at  1.33  feet.  Instead,  the  transition 
occurred  at  0.83  feet.  This  corresponds  to  a  free  stream  turbulence 
Intensity  of  about  0.4  percent.  The  corona  was  turned  on  leaving  the 
premature  transition  unexplained  for  the  moment. 

Two  observations  can  be  made  from  the  test  results  shown  on 

Figure  35: 

(1)  The  location  of  transition  was  unaffected 

by  the  corona  discharge,  both  with  and  without 
the  controlled  frequency  component  (1215  cps). 

(2)  A  decrease  in  impact  pressure  behind  the  last 
wire  was  present  when  the  discharge  was  on, 
both  with  and  without  the  controlled  frequency 
component . 

In  an  effort  to  determine  the  cause  of  the  premature  transi¬ 
tion,  the  wires  were  removed  and  the  flow  adjustment  doors  set  fully 
open  to  lower  the  flow  speed.  If  the  transition  point  were  now 
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located  downstream  of  1.33  feet,  the  premature  transition  could  be 
attributed  to  presence  of  the  wires.  If,  however,  the  transition 
point  was  between  0.83  feet  and  1.33  feet,  the  wires  were  probably  not 
the  cause.  A  test  was  conducted  at  153  fps  and  transition  was  located 
at  1.17  feet.  Three  wires  were  immediately  installed,  at  distances  of 
1.063,  1.112,  and  1.16  feet  from  the  leading  edge,  and  a  test  was  con¬ 
ducted  at  153  fps  with  and  without  current.  Provision  was  not  made  to 
change  the  wire  spacing  but  the  excitation  frequency  was  Increased  to 
compensate  for  the  change  in  location  on  the  stability  diagram. 

Figure  4.  These  observations  may  be  made  from  the  test  results  shown 
in  Figure  36: 

(1)  Transition  was  not  tripped  by  the  wire. 

(2)  The  location  of  transition  was  unaffected  by 
the  corona  discharge,  both  with  and  without 
the  controlled  frequency  component 

(3)  The  impact  pressure  dip  marking  transition 
was  accentuated  in  the  presence  of  the  dis¬ 
charge,  both  with  and  without  the  controlled 
frequency  component. 

At  the  conclusion  of  this  series  of  tests,  the  reason 
explaining  the  premature  transition  at  175  fps  was  made  evident.  A 
guard  rail  had  been  recently  installed  just  outside  the  entrance  to 
the  tunnel.  The  rail  may  be  seen  in  Figure  16.  The  obvious  effect 
of  the  rail  was  to  increase  the  free  stream  turbulence.  This  caused 
the  normal  transition  to  move  from  1.33  feet  to  0.83  feet.  The  rail 
was  renoved  for  the  final  series  of  tests  now  to  be  described. 
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Following  the  experimental  procedure  outlined  above,  a 
final  scries  of  corona  tests  was  conducted  at  174  fps.  The  first  test 
was  performed  with  no  wires  installed.  The  static  pressure  distribu¬ 
tion  for  this  test  is  shown  in  Figure  26.  Transition  was  located  at 
1.25  feet.  Four  wires  (0.008-inch  diameter)  were  then  installed  from 
1.063  to  1.209  feet  from  the  loading  edge  at  a  spacing  of  0.048  feet. 
Another  test  was  conducted  at  174  fps  with  no  current  on.  The  transi¬ 
tion  was  essentially  unchanged,  1.27  feet.  For  the  next  part  of  the 
test,  the  electrical  disturbance  was  activated.  The  corona  parameters 
were  as  follows: 

Current,  1  ma 

Voltage,  9  kvdc  plus  2  kvac  peak-to-peak 

Frequency,  1215  cps. 

The  transition  point  shifted  downstream  to  1.42  feet.  At  1.27  feet, 
the  pressure  readings  were  now  steady  and  transition  was  not  apparent. 

To  confirm  the  transition  shift,  the  current  was  later  turned 
completely  off.  Impact  pressure  readings  were  taken  immediately  at 
1.42  feet  and  1.27  feet.  Fluctuating  pressure  readings  at  a  given 
distance  from  the  leading  edge  signify  the  presence  of  transition,  as 
discussed  above  under  DETECTING  TRANSITION.  Once  again,  the  pressure 
readings  were  steady  at  1.42  feet  and  fluctuating  at  1.27  feet,  indi¬ 
cating  that  the  transition  had  jumped  back  upstream.  The  data  for 
the  tests  Just  described  are  presented  in  Figure  37.  Leading  edge 
data  have  been  omitted  to  preserve  scale.  All  sets  of  data  terminate 
at  33.2  feet  of  water  at  the  leading  edge. 
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During  the  test,  three  sweeps  of  excitation  frequency  were 
made.  The  sweeps  were  made  at  probe  positions  of  2  feet,  1.42  feet, 
and  1  foot  from  the  leading  edge.  At  each  frequency,  shown  on 
Figure  38,  one-half  minute  was  waited  before  recording  the  impact 
pressure.  Shifting  between  readout  frequencies  was  accomplished  by 
slowly  adjusting  the  oscillator  with  the  vernier  knob.  In  this 
manne*-,  it  was  felt  that  any  resonant  or  unusual  effects  would  be 
revealed.  Such  was  not  the  case,  for  the  transition  apparently  re¬ 
mained  at  1.42  feet  regardless  of  the  frequency.  This  was  also  true 
for  zero  frequency,  or  straight  d-c  corona. 


I 
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DISCUSSION  OF  RESULTS 


The  effect  of  electrical  disturbance  on  transition  will  be 
discussed.  In  the  preceding  section,  it  was  demonstrated  that  corona 
discharge  may  be  used  to  forestall  transition  under  certain  conditions. 
Several  questions  arise  regarding  the  phenomenon.  What  is  the  mecha¬ 
nism  by  which  transition  was  forestalled?  Why  was  the  transition 
insensitive  to  tht  electrical  disturbance  frequency?  Lastly,  why  was 
the  transition  affected  at  the  flow  speed  of  175  fps  and  unaffected  at 
153  f ps? 

In  an  effort  to  disclose  the  mechanism  influencing  transition, 
several  characteristics  of  the  electrical  discharge  will  be  brought 
forth.  First,  for  a  constant  sustained  discharge,  the  total  velocity 
of  the  ions  is  the  sum  of  air  velocity  and  the  velocity  of  the  ions  in 
the  air. 


v 


t 


(46) 


The  velocity  of  the  ions  in  air  is  given  by  the  product  of  the  field 
strength  and  the  ion  mobility 


ion 


K  E 


(47) 


For  a  positive  wire  corona,  it  may  be  assumed  that  only  positive  ions 

exist  in  the  major  portion  of  the  wire-to-plate  gap.  For  air  at  normal 

2 

pressure  and  temperature,  Reference  16  gives  a  value  of  0.00146  ft  / 
volt-sec  for  the  mobility  of  singly  charged  positive  ions.  If  the 
field  between  the  parallel  wires  and  the  plate  is  approximated  as  that 
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between  two  parallel  plates,  the  field  voltage  drop  may  be  taken  as 
half  the  wire-to-plate  voltage  because  of  the  sharp  voltage  peak  near 
the  wires.  Therefore,  the  field  strength  is  approximated  by 

E  .  I10,000/2)^vol_t:8x_12  tea  »  107,000  volts/ft  .  (48) 

Using  Equation  47  gives  an  ion  velocity  of  156  fps. 

The  ion  velocity  is  close  to  the  air  flow  velocity  outside 
the  boundary  layer.  The  ions  are,  therefore,  impinging  on  the  boundary 
layer  at  approximately  a  45-degree  angle.  The  relative  absolute  veloc¬ 
ity  of  the  neutral  air  molecules  in  the  corona  wind  will  be  about 
(21) 

5  fps  .  Thus,  the  neutral  air  molecules  impinge  on  the  boundary 
layer  at  a  shallow  angle.  The  static  pressure  Increase  at  the  plate 
that  would  be  associated  with  a  direct  corona  wind,  that  is  if  the 
tunnel  flow  velocity  were  zero,  is  of  the  order  of  only  0.001  feet  of 
water,  according  to  actual  experiments  performed  in  Reference  12. 
However,  it  demonstrated  that  such  a  corona  wind  is  capable  of  greatly 
depressing  a  convective  boundary  layer  with  associated  drastic  in¬ 
creases  in  local  heat-transfer  coefficient.  It  is  postulated  that 
similar  boundary-layer  depression  is  occurring  in  the  present  experi¬ 
ment  to  a  lesser  extent.  In  the  absence  of  a  boundary -layer  mouse, 
this  could  not  be  verified. 

The  boundary-layer  depression  manifests  itself  by  causing 
high  velocity  streamlines  to  be  brought  closer  to  the  surface  of  the 
plate  in  the  region  adjacent  to,  and  downstream  of,  the  wires.  This 
may  be  seen  in  Figure  36  independent  of  transitional  effects. 
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Hypothetically,  if  the  wires  are  positioned  immediately  upstream  from 
the  transition  point,  boundary-layer  depression  would  imply  two  things. 
First,  the  boundary  layer  is  being  physically  confined  to  the  plate. 
Second,  the  thinning  boundary  layer  will  have  associated  with  it  a 
local  negative  static  pressure  gradient  in  the  flow  direction.  It  is 
well  known  that  transition  is  retarded  by  a  favorable  pressure  gradient. 
Thus,  boundary -layer  depression  may  be  the  mechanism  explaining  the 
effect  of  the  electrical  disturbance  on  transition. 

The  lack  of  response  to  the  controlled  frequency  component 
of  the  electrical  disturbance  may  be  explained  in  several  ways.  It  is 
felt  that  one  of  the  primary  reasons  is  Insufficient  intensity  of  the 
a-c  component.  Current  contributions  associated  with  the  pulsating 
disturbance  were  an  average  of  20  ^a.  This  is  2  percent  of  the  total 
corona  current.  Another  reason  may  be  Improper  matching  and  phasing 
of  the  disturbance  with  local  predominant  boundary -layer  oscillation. 
Hot-wire  equipment  would  be  required  to  determine  the  actual  wave 
properties  of  laminar  boundary-layer  oscillations  present  in  the  flow. 
Without  this,  choosing  a  disturbance  from  the  theoretical  neutral 
stability  curve  could  never  precisely  match  a  predominant  oscillation. 
Furthermore,  It  Is  not  known  by  what  mechanisms  the  amplification 
process  would  be  disrupted  even  if  proper  matching  and  phasing  could 
be  accomplished.  However,  it  can  be  concluded  from  the  experiment 
that  the  controlled  frequency  component  of  the  disturbance  did  not 
have  a  detrimental  effect  on  the  amplification  process.  This  would 
cause  premature  transition. 
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In  Test  17,  transition  was  forestalled  by  the  action  of  the 
corona  while  in  Test  IS  it  was  not.  In  both  tests,  the  corona  wires 
were  positioned  immediately  upstream  of  the  normal  transition  location. 
The  explanation  lies  in  the  difference  of  free  stream  turbulence  in¬ 
tensity  between  the  two  series  of  tests.  Recall  that  Tests  14  and  * 
were  conducted  in  the  presence  of  the  outside  guard  rail  while  Test  17 
was  not.  As  mentioned  above,  the  deduced  free  stream  turbulence  in¬ 
tensity  associated  with  Test  IS  was  0.4  percent  and  that  for  all  other 
tests  was  approximately  0.3.  This  30  percent  Increase  in  background 
turbulence  might  well  be  the  swaying  factor  in  determining  whether  or 
not  the  corona  discharge  is  capable  of  influencing  transition.  From 
this  it  may  be  suggested  that  future  investigations  be  conducted  at 
the  lowest  possible  levels  of  free  stream  turbulence. 
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SUMMARY  OF  CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  resultant  from  this  experimental  investiga¬ 
tion  are  summarized  below  in  chronological  order. 

(1)  To  assure  proper  development  of  a  laminar  boundary 
layer,  on  a  flat  plate  with  an  asymmetric  sharp  nose,  requires  the 
use  of  flow  "blockage"  techniques  on  the  active  side  of  the  plate. 

(2)  A  near-surface  impact  probe  has  proved  to  be  a  reliable 
and  convenient  tool  for  detecting  boundary -layer  transition  providing 
the  probe  tip  size  is  carefully  chosen.  It  was  found  that  a  tip 
diameter  approximately  equal  to  1/3  the  laminar  boundary-layer  thick¬ 
ness  at  transition  (5fcr)  was  quite  satisfactory  from  the  aspects  of 
both  response  and  sensitivity. 

(3)  It  was  found  that  rows  of  wires  oriented  at  right 
angles  to  the  flow;  suspended  at  a  distance  of  9  x  6fcr  from  the  flow 
surface  and  extending  from  the  normally  laminar  region  deep  into  the 
turbulent  region,  cause  transition  to  occur  at  the  wire  farthest  up¬ 
stream.  Wires  similarly  suspended  in  only  the  pre-transition  region 
do  not  trip  transition. 

(4)  Transition  is  sensitive  to  free  stream  turbulence.  Small 
flow  obstructions  outside  the  intake  zone  of  the  wind  tunnel  may  cause 
premature  transition. 

(5)  Under  conditions  of  low  free  stream  turbulence,  a  steady 
corona  discharge  applied  to  the  boundary  layer  in  the  pre-transition 
region  can  forestall  transition.  Boundary-layer  depression  is  the 
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likely  mechanism  involved.  This  was  accomplished  at  a  mean  flow  speed 
of  175  fps. 

(6)  Superposition  of  a  controlled  frequency  component  on 
the  steady  corona  discharge  did  not  affect  the  result  stated  in  (5), 
nor  did  it  promote  premature  transition.  A  wide  range  of  frequency 
was  scanned.  The  failure  of  the  controlled  frequency  excitation  to 
interact  with  the  selective  boundary-layer  amplification  process  is 
attributed  to  lack  of  intensity  and  precise  matching  of  wave 
parameters. 

In  the  light  of  the  experimental  findings,  it  appears  logi¬ 
cal  to  continue  further  investigation  of  electric  effects  on  boundary- 
layer  transitions.  Some  recommendations  regarding  the  course  and 
methods  of  such  experimentation  are  presented. 

(1)  The  free  stream  turbulence  Intensity  should  be  low, 
certainly  below  0.3  percent. 

(2)  A  hot-wire  anemometer  should  be  employed  to  probe  the 
boundary  layer  for  thickness,  laminar  oscillations,  and  transition. 

(3)  In  addition  to  detecting  the  natural  disturbances 
present  in  the  boundary  layer,  provision  should  be  made  to  determine 
the  effects  of  the  forced  perturbations.  Local  pressure  changes 
caused  by  the  action  of  the  corona  may  be  sensed  with  suitably  placed, 
sensitive  instrumentation. 

(4)  Means  for  attaining  a  controlled  frequency  disturbance 
of  greater  intensity  should  be  sought. 
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(5)  With  respect  to  the  present  electrical  configuration, 
provision  to  vsry  the  wire  spacing  and  positions  would  be  most 
desirable.  The  effects  of  wires  positioned  near  the  leading  edge 
and  into  the  transition  region  should  be  explored. 

(6)  Different  electrode  configurations  should  be 
investigated. 

(7)  Slower  flow  speeds  may  be  beneficial  in  highlighting 
sensitive  ion-fluid  interactions. 
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Reynold*  Number ,  R  g* 

FIGURE  3.  NEUTRAL  STABILITY  CURVE  IN  TERMS  OF  REYNOLDS  NUMBER  AND 
VELOCITY  PARAMETER 
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Reynolds  Number,  Re  x 

FIGURE  4.  NEUTRAL  STABILITY  CURVE  IN  TERMS  OF  REYNOLDS  NUMBER  (Re„) 
AND  FREQUENCY 


59 


FIGURE  7.  SECTION  VIEW  OF  FLAT  PLATE  LEADING  EDGES 


FIGURE  8.  FLAT  PLATE  STATIC  PRESSURE  TAPS 
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FIGURE  10.  DOWNSTREAM  VIEW  OF  TEST  SECTION 


FIGURE  12.  ELEVATION  VIEW  OF  WIND  TUNNEL 


Ratio  Of  Local  Dynamic  Pressure  To  Full  Dynamic  R'essure 


FIGURE  14.  HORIZONTAL  DISTRIBUTION  OF  DYNAMIC  PRESSURE 
AT  CENTERLINE  OF  TEST  SECTION 


Distance  From  Floor,  feet 
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FIGURE  15.  WALL  BOUNDARY  LAYER  AT  THE  VERTICAL 

CENTERLINE  OF  WIN!'  TUNNEL  TEST  SECTION 
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FIGURE  16.  WIND  TUKKEL  AIR  INTAKE  AREA 
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FIGURE  17.  TURBULENCE  SCREENS  HEAR  WILT)  TUIii.'EL  AIR  INTAKE  AREA 
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FIGURE  18  IMPACT  PRESSURE  IN  THE  BOUNDARY  LAYER  AT  A  SMALL 
CONSTANT  HEIGHT  ABOVE  THE  SURFACE 
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FIGURE  22  PULSATING  POTENTIAL  IMPOSED  ON  CORONA  WIRES 


75 


i 

I 

I 

I 

I 

I 


I 


FIGURE  23.  CIRCUIT  DIAGRAM 
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FIGURE  2k.  ELECTRICAL  EQUIPMENT  AID  TEST  SECTION  EXTERIOR 
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FIGURE  25.  MANOMETER  BOARD  FOR  INDICATING  STATIC  PRESSURE  DISTRIBUTION 
OVER  PLATE 
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FIGURE  26.  TYPICAL  STATIC  PRESSURE  DISTRIBUTION 
OVER  PLATE  SURFACE 


FIGURE  27.  TUNNEL  CONTROL  CONSOLE  SHOWING  BAROMETER  AND 
FLOW  CONTROL  MANOMETERS 
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FIGURE  30.  NEAR  SURFACE  IMPACT  PRESSURE  WITH  SHARP 


LEADING  EDGE 
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LAYER  TRANSITION 


FIGURE  33.  FREQUENCY  PARAMETER  VERSUS  REYNOLDS  NUMBER  FOR  §  *0.24 
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Neor-Surfoce  Impact  Pressure,  feet  of  water 
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FIGURE  35.  CORONA  TEST,  FLOW  SPEED  175  FPS 
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FIGURE  36  CORONA  TEST,  FLOW  SPEEO  153  FPS 
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FIGURE  38.  IMPACT  PRESSURE  AS  A  FUNCTION  OF  THE  FREQUENCY 
OF  THE  ELECTRICAL  DISTURBANCE 
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